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Implications of the internal symmetries on the dynamics of the trapped two-component atomic 
vapors are discussed. In the cases of *^Rb (bosons) as well as of ''"K (fermions) trapped in the 
I two hyperfine states, the intrinsic su{2) symmetry can be realized with a very good precision. 

. Such a symmetry protects the global operators, which are the generators of the symmetry, from 

any decoherence. The case of boson-fermion mixture is discussed as well. The role of external 
factors, breaking the symmetry, in inducing the decoherence of the global operators is considered. 
It is shown that, the loss of the correlations is not faster than the rates of the induced heating or 
losses, provided the noise is characterized by short correlation length. The case of the extrinsic 
long-ranged fluctuations is also considered. Intrinsic mechanisms of decoherence of the correlators 
of the condensate operators of the two-component condensate are analyzed. The atomic echo is 
discussed as a test for the reversibility of the phase diffusion effect. The intensity and the profile of 
. the multiple echo are calculated numerically as well as analytically. 
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I. INTRODUCTION 



Successful achievements of confinement and cooling to the state of degeneracy of Bose and Fermi |^ atomic gases 
'■^ , have created unprecedented variety of possibilities for studies of the many-body effects. The fundamental issue of 

■ emergence of irreversibility can now fully be addressed experimentally and theoretically in the wide range of situations. 
O [ Studies of the dissipation of the normal modes of the confined Bose-Einstein condensates (BEC) Q began immedi- 

I ately after BEC has been realized. Very recently, the data on the damping in the Fermi mixtures has been reported 

■ as well g. 

The Ramsey spectroscopy study of the temporal correlations of the relative phase of the two-component ®'^Rb BEC 
K*' has been conducted by the JILA group j|] . No decay of the global relative phase has been detected on the time scale 
, of the experiment < 100ms. The question was then posed in Q as to why the phase correlations are so robust despite 
' an apparent fast relaxation of other degrees of freedom such as, e.g., the relative motion of the two condensates. 
I Many intriguing questions are associated with the finiteness of the number of bosons N forming BEC in the trap. 

■ It has been realized that the global phase may experience collapses and revivals as long as N is finite Q . 

T-H Symmetries play a very special role in the many-body dynamics. As was discovered by Pitaevskii and Rosch a 

i 2D gas trapped by the oscillator potential and characterized by the contact two-body interaction can be described by 
' the dynamical symmetry 5*0(2, 1), and this results in the non-decaying breathing mode. The symmetry protects the 
^ mode from the decay. In our recent work it has been found that, in the case of the ^^Rb mixture characterized 
^ by the proximity to the intrinsic su{2) symmetry the dynamics of the global generators is also protected by 
' , this symmetry from the intrinsic decoherence. In the case of the broken symmetry, the rate of the decoherence is 
(— ( controlled by the deviations from the su{2), and, accordingly, the decoherence is slow as long as the deviations are 
Q small. This answers the question Q about the robustness of the phase dynamics. 

O In another our work the nature of relaxation of the condensate operators was discussed. I has been shown 

that, if the number of bosons A'o in the 3D equilibrium condensate is macroscopic, the decoherence of the condensate 
time-correlators is rather a reversible dephasing than a true irreversible decay. In some sense, the Bose permutation 
symmetry, representing the indistinguishability of particles which populate the lowest (condensate) single particle 
state, protects the condensate operators from the irreversible loss of memory. It was suggested to test this by the 
echo-type experiment. Thus, even though the intrinsic su{2) symmetry can be broken, the dynamics of the condensate 
operators remains reversible. 

In the present work, we, first, extend our analysis |^ to the case of the two-component ferniion mixture [Q, and 
will show that, if the trapping potentials for the components are identical, the global quantities - the generators 
of the symmetry formed by the Fermi operators - do not decay in spite of the presence of the interaction between 
the components. We also analyze the case of the boson-fermion mixture, which is close to the su(2) symmetry 
with respect to the mutual (formal) transformation of bosons and fermions into each other. In this situation, the 
decoherence of certain global operators, which, however, are not the generators, is suppressed as well. 

Second, we analyze the problem of the decoherence in the case of deviations from the su(2) symmetry in the 
two-component BEC. Extrinsic and intrinsic mechanisms of the decoherence are considered. For the extrinsic factors. 
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we concentrate on the interaction with the hot background gas, which induces heating and losses (or deposition) . The 
heating is modeled by introducing a random white noise potential, which is characterized by some spatial correlation 
length Lp. We show that the ratio of the phase decay rate t~1 to the heating rate r^^ is essentially given by the 
factor Lp/L^, where Lc stands for the thermal length of the trapped atoms at the BEC transition temperature Tc- 
Obviously, Lp is given by the thermal length of the background gas, which is <C Lc- This implies that the phase 
memory loss turns out to be not faster than the heating rate. Similarly, incoherent losses from the trap are unable to 
produce the phase decay faster than the rate of the losses. 

We also analyze the role of the instrumental effects, which may produce the long range noise breaking the intrinsic 
symmetry. In this situation, the phase decay rate can be significantly enhanced, if compared with the corresponding 
rate of disrupting the thermal equilibrium within the atomic cloud. However, we estimate that the phase decay rate 
is still too small in the current traps, provided the instrumental effects are on the level of the equilibrium room 
temperature thermal noise. 

Regarding the intrinsic factors, we extend our analysis to the case of the two-component BEC, and show that 
the intrinsic factors do not induce the irreversible decay of the correlators of the condensate operators. At finite 
temperatures T, the phase exhibits reversible phase diffusion effect similar to that analyzed previously at T = by 
many authors In fact, in any exact many-body eigenstate, the evolution of the condensate operators is given by 
the relative chemical potential. Averaging over the thermal ensemble produces dephasing, whose rate is determined 
by the fluctuations of the chemical potential. At T = 0, these fluctuations are produced by the averaging over the 
initial coherent state and/or shot noise, and, at T ^ 0, the thermal averaging induces an additional contribution. 
Thus, at finite T the nature of the global relative phase diffusion remains, practically, the same as that at T = 0, 
provided intrinsic factors only contribute to the decoherence. 

Finally, we analyze the atomic echo effect, suggested previously for testing the reversibility of the phase diffusion 
pO|,p[, in the context of the JILA experiment g. The strength of the echo is evaluated. When the dephasing is 
dominated by the shot noise in the total number of bosons in the trap, the echo can be as large as 100%, for the 
TT-pulse employed for the time reversal When the phase diffusion is dominated by the self- interaction effects at 
T = 0, the echo about 50% can be achieved by a weak time-reversal pulse, which transfers coherently a (relatively) 
small number of bosons between the components. An analysis of the echo at finite temperature is presented as 
well. Our conclusion is that the echo, although suppressed, can still be observed at finite temperatures, provided the 
population of the condensate state is dominant. 



II. INTRINSIC SYMMETRIES AND THE DYNAMICS OF THEIR GENERATORS IN THE 

TWO-COMPONENT ATOMIC MIXTURES 

Here we consider simplest cases of the two-component mixtures: bosons of sort 1 + bosons of sort 2; fermions of 
sort 1 -|- fermions of sort 2; and one component bosons + one component fermions. 

A. Two-component Bose gas 

First, for the purpose of consistency, we will briefly outline the results Q for the two-component Bose mixture. We 
employ the following two-component Hamiltonian 

H = J dx{*I(i/i-|)*i + M/t(i72 + |)vI/2 + (^^vI,tvi,^+i/.c) + 

(1) 

+ ^^[ai'I'l^'J*!*! + a2^'J*J*2*2 + 2ai2^'I*J^'2*i]}, 

where ^1*^ is the second quantized Bose field of the i-the component {i = 1,2); Hj — —h^V'^/mj + Uj{x.) stands 
for the one-particle part which includes the kinetic energy and the trapping potential; the quantity eg = const is the 
detuning of the external field h^{t) (taken in the rotating wave approximation ) from the resonance between the 
components; Q(t) stands for the corresponding Rabi frequency treated as an envelope of the rf-pulses; The binary 
collision terms in (|l]) are taken in the contact form, with ai, 02, ai2 being the corresponding scattering lengths. 

The intrinsic su{2) symmetry corresponds to the situation Hi = H2, ai — a2 = ai2. The generators of this 
symmetry 
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Iz = lJ dxC^^^-a - ^-I^i); I+= J dx^t^i; = J dyL^^2 (2) 

represent the su{2) algebra of the angular momentum operators. Indeed, it is easy to see that these operators obey 
the standard su{2) commutation relations 

[4,/+] = /+, = [/+,/_] = 24, (3) 

provided the field operators obey the Bose commutation rule [^'i(x), 5'j(x')] — (5y(5(x — x'). The operator Iz = 
{N2 — Ni)/2, with standing for the total number of bosons in the first and the second components, respectively, 

can be viewed as the z-component of the angular momentum operator; = (/+ +/-)/2, ly = (/+ — /_)/2i arc the 
x,y-components, respectively. 

If the condensate is present, the operators I± carry the information about the relative global phase (p (we may call 
it simply "phase" ). Indeed, in this case, ^'j w y^^jW '^^'^^ 7 where (pj is the global phase of the jth component, and 
V stands for the effective volume occupied by the condensate. Thus, /+ « y/ThT^e^'^ , cp = ipi — (p2- Accordingly, 
the correlator (/+(i)/_(0)), where the mean is taken with respect to some set of initial states, carries the information 
about the time correlation properties of the phase. 

In the case of the exact intrinsic symmetry, the generators (H) obey linear Heisenberg equations 

= thi+ = -eoi+ + nn{t)iz, tnL = e„i--nn*{t)iz, (4) 

where the su{2) commutation relations (^) have been employed If n{t) and eg are not subjected to any noise, 
the solution for the correlator (/+(t)7_(0)) exhibits no decoherence despite the presence of the two-body interaction 
in the Hamiltonian (|l|). It should also be noted that Eq.(^ are exact, and they do not require a presence of the BEC. 

The decoherence will arise, if cither the one-particle parts Hi 2 are different or the scattering lengths deviate from 
the condition oi = 02 = ai2. The nature of the decoherence in the case of the BEC will be discussed later. 



B. Two-component Fermi mixture 



In the case of the two-component Fermi gas, the su{2) intrinsic symmetry can hold as well. For example, trapping 
''"K in the optical trap results in the identical potentials for both components This, in fact, implies that the 
intrinsic su{2) symmetry holds for any value of the inter-component scattering length. Indeed, in the Hamiltonian 
, the Bose operators \l/i,2 are to be replaced by the Fermi operators Fi^2 obeying the anti-commutation rule 

{F,(x),7;t(x')} = 5,;,J(x-x'), (5) 

where {...,...} denotes the anti-commutator. Accordingly, the terms ~ ai, 02 vanish in Eq.(0). Introducing the 
su{2) generators (|^), where the Bose-operators are replaced by the Fermi-operators, the su{2) commutation relations 
follow. Finally, the equations of motion for the generators turn out to be Eq.lQ). Thus, in the case of the two- 
component fermion gas, the dynamics of the generators is always linear regardless of the value of the scattering length 
ai2, as long as the one-particle Hamiltonians Hi 2 are the same. We should, however, note that the higher order 
scattering harmonics change this situation, and the corresponding restrictions should be imposed on the scattering 
amplitudes of the higher orders (p-,f- waves etc). These processes, however, are insignificant at low temperatures and 
at low densities. 



C. Bose- Fermi mixture 



In the case of the mixture of, e.g., one-component Bose and Fermi gases, the respective fields ^'(x), F(x) are 
introduced, and the Hamiltonian takes the form 

dx{'i>\Hi - + fUH2 + + ^^[ai^t^t^* + 2ai2'i'''F^Fm, (6) 
2 2 m 

where it was taken into account that no interconversion between bosons and fermions is possible, and that the 
self-interaction between the one-component fermions vanishes in the s-wave approximation. 



3 



Let us analyze a special situation - Hi = H2, ai = 012. This corresponds to the (accidental) intrin- 
sic su{2) symmetry between the bosons and the fermions. Indeed, the interaction term can be rewritten as 
a^(\I/t\i/)2 _|_ 2ai2^'''^F'' F = ai[*^\l' + F'' F]^ in this case, because [F^ F]'^ vanishes (apart from a trivial term 
~ F^F, which simply renormalizes the constant eq in Eq.(|6|)). Now, if Hi = H2, the formal su{2) transformation 
acting on the "spinor" made of F retains the Hamiltonian intact. 

We introduce the operators 

Iz^^ J <ix(*t^ - FtF), I+= J dx-^^F, I-"^ J dxFt^-, 7 = ^ dx{¥-^ + F^F). (7) 

It should be noted that these operators do not form a closed algebra, as can easily be verified by commuting them 
with each other. Instead, they obey the following relations 

= /+, = -/_, {/+,/-} - /, [/.,/] = [/±, /] - 0. (8) 

It can be verified that the Heisenberg equations of motion for the above operators are also linear. Specifically, 

t =0, / = 0, ini±^TenI±, (9) 

as long as the above condition of the intrinsic su{2) symmetry holds. Thus, the situation of the Bose-Fermi mixtures 
can also exhibit the dissipationless dynamics of certain global quantities (which are not, however, the generators of 
the symmetry) . It should be mentioned that any deviations from the intrinsic symmetry will result in the decoherence 
of these global quantities, and the rate of such a decoherence will also be controlled by the proximity to the symmetry. 

III. THE NATURE OF DECOHERENCE IN THE CASE OF THE BROKEN SU{2) SYMMETRY IN THE 

TWO-COMPONENT BOSE-MIXTURE 

Here we will discuss various factors which break the intrinsic su{2), and thereby introduce the decoherence in, e.g., 
the correlator 

p,,(x,x',i,<') = (*I(x,t)*,(x',t')), t>t'. (10) 

Hereafter we consider the Bose statistics only. 

First, we address the decay induced by extrinsic factors such as: 1) collisions with hot background gas, which 
result in heating of the confined cloud and in losses from the trap; 2) Thermal noise of the trapping potential, which 
corresponds to temperatures much higher than the confined gas. 

We raise a general question: Under what conditions can the loss of the memory, produced by the extrinsic factors, 
be much faster than the corresponding time of disrupting the equilibrium in BEG! In other words, if is the time 
of the heating induced by some extrinsic factor, and Tirr stands for the time of the irreversible loss of the time 
correlations in (|l^) due to the same factor, can it be that the condition 

Tirr < Th (11) 

holds? In the same sense, if tl stands for a typical life-time of the confined cloud, which is being subjected to losses, 
can the condition ([III ) hold, where the role of Th is replaced by tl ? 

Below we will show that the range of the spatial correlations of the extrinsic factors plays a crucial role in fulfilling 
the condition (|ll|). In fact, the case 2) turns out to be the most efiicient in destroying the time correlations without 
introducing significant disturbances into the system. 

A. Rates of disrupting the equilibrium, and the loss of the phase memory 

The purpose of the following is to outline the framework within which the effects of environment on the confined 
cloud can be analyzed, and to obtain general criterion allowing to compare the rates of disrupting the equilibrium 
and of the decoherence. In this section, we will consider the effects of interaction with the background (hot) ga^. In 



^ We believe that the conclusion ensuing from this analysis is general 
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order to do this consistently, the corresponding interaction term should be added to the Hamiltonian (|]). We choose 
it in the form 



Hr 



ij,kl 



J dxg,,-fci*^4Si*„ (12) 



where the summation runs over all the components of the confined cloud and of the background gas described by the 
fields Bi] QijM are the corresponding interaction constants. 

The term (|lj) describes two basic effects [|2| : a) a fast particle strikes a confined boson and ejects it from the trap; b) 
the background particle does not eject any boson from the cloud and, instead, imparts a substantial energy producing 
heating in the trap, and, possibly, exchange between the components. The term (O) can be employed in order to derive 
the equation for the one-particle density matrix (OPDM) /9y (x,x',t) = (4'|(x, t)^'j(x', i)) = (x, x', t, i) as well 
as for the correlator ( p^ ) . The OPDM equation becomes a set of the kinetic equations in the Wigner representation. 
Here we will not implement this to full extent. Instead, we will analyze an exactly solvable model which captures 
basic features of the effect b), and then we will present results for the effect a). 

The field = gij^kiB\Bi in Eq.(|l2|) can be viewed as some random mean field potential acting on the confined 
cloud. As long as the background gas is hot, this effective random potential can be considered as a white noise. As it 
will be seen later, this assumption allows obtaining the OPDM equation exactly in the case when collisions between 
the confined bosons are ignored. 

In order to simplify further analysis, we ignore the exchange of the components induced by the collisions^, and 
consider just a one-component random potential ^ which satisfies the condition 

(e(x,i)e(x',0>=^(i-W|x-x'|), (13) 

where C(|x — x'|) stands for some smooth function which decays to zero on some typical length Lp. This function 
describes spatial correlations of the noise. 

It should be noted that the above model is not suitable for treating a thermal equilibrium between the background 
and the cloud. Indeed, the external time dependent potential may produce unlimited heating, similarly to the case 
of the Brownian motion, if no friction term counterbalancing the external noise is introduced. In order to treat the 
equilibrium, a fully Hamiltonian analysis based on Eq.(|l^) should be implemented. This, however, is not required 
because the confined cloud is essentially out of the equilibrium with the background, and we are estimating the rate 
of disrupting the equilibrium in the cloud. 

The Hamiltonian is taken in the form (|l]), where the noisy potential should be added. We choose ai = 02 = ai2 = 0. 
Furthermore, for the sake of simplicity, we consider no trapping potential and set without the loss of generality 
£0 = 0, = , so that the Hamiltonian takes a form 

F = y" dx{^'ti/„^,^ + ^t^Q*2+ei«'l*i +6*^*2]}, ei--e2 = ^, (14) 

where Hq = — |^V^, and we have included the interaction with the random potential ^ through the term breaking 
the intrinsic su{2) symmetry, because, according to the preceding analysis, no decoherence of the global operators 
can be induced in the case of the intrinsic symmetry ( ^1 = ^2 hi the considered situation). 



In the presence of the noise, the meaning of the averaging in Eq.(lO) should be specified. Let us assume that some 
initial state |i = 0) was created at the time moment t — Q, and the noise affects the following evolution at later times 
i > 0. Thus, the averaging should be performed, first, over the initial state (or a set of states with some weight), and, 
then, the averaging over the gaussian noise is to be done. 

The OPDM characterizes the rate of the heating. Indeed, the mean kinetic energy of the particles is defined as 

/ dxi/ox'Pi,:(x,x',i), x' ^ x, (15) 

i=l,2 

where, the coordinate dependence -ffox' hr the kinetic energy operator Hq indicates that this operator acts on the 
coordinate x', and, then, one should set x' = x and integrate over x. 



while complicating the calculations, inclusion of the random exchange does not change the main conclusion 
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Equations for the OPDM follow from the Heisenberg equations ih'^a = {Hq + ^a)^a- Employing the Furutzu- 
Novikov theorem (see details in Appendix A), we obtain the exact equations for the OPDM as 



2i 

ihpu{x,x.' ,t) = [-i/ox + Ho^']pu{x,x.' ,t) - -^[C(O) - ^dx - x'|)]pii(x, x', t); 



(16) 



ihpi2{^, x', t) = [-iJox + Fox']pi2(x, x', t) - -r [C(0) + C(|x - x'|)]pi2(x, x', t). 



(17) 



It should be noted that reinstating the inter-particle interaction will result in a chain of coupled equations for the 
correlators, where the mean involves averaging over the noise. On the each next step, the decoupling of the noise term 
can be done by means of employing the Furutzu- Novikov theorem p3| . This will lead to finding the corrections to the 
collision integral due to the noise. Practically, it is enough to neglect such corrections. Accordingly, Eqs.([l6|), ( |T^ ) 
can be modified by simply adding the collision integrals obtained with no noise present. Inclusion of such integrals, 
however, does not affect the main result. 

We choose the initial state as a uniform condensate characterized by finite correlations between the components. 
Accordingly pij(x, x',t = 0) = p^^"* = const for all the components. As time proceeds, this state becomes destroyed. 
Indeed, the solution of Eqs.([l^), corresponding to the chosen initial condition is 



Pii{x,x',t) = e 



^2[C(0)-C(|x-x'|)]t/n^ (0) 



pi2(x,x',i) = e-2[c^(o)+c^(|x-x'l)]tA^^(o). 
The increase of the kinetic energy, then, follows from Eq. (jlj) as 

if = X(<) = --VLoC(|x|)t, Kit)^^^r 



(18) 
(19) 
(20) 



(21) 



where the relations pulx, x, t) — p^^\ {p^-^]^ + ^220^ = ^ i with V being the total volume of the system, have been 
employed. We have also defined the noise correlation length Lp through the relation V^^oC'dxl) = — C(0)/Lp, and 
imposed the condition Vx=oC'(|x|) ~ 0. 

We note that Eq.(|2^) accounts for the loss of the correlations between the components with the effective rate 
« 2(C(0) + C(x — x' ))/n^. This rate can be obtained more accurately from the equation for the correlator 
Pi2(x, x', t, t' = 0), which describes the decay of these correlations and carries the information about the relative phase 
(below Tc). Applying the Furutzu-Novikov theorem jl^ again, and assuming that the initial state is independent of 
the noise, we find 



pi2(x,x',t,i' = 0) = e 



-c(o)tA^^(o)^ i>0. 



(22) 



Thus, the rate of the loss of the phase memory is 



C(0) 



irr ^2 



(23) 



which is a factor of 2-4 different from the estimate obtained from the OPDM. 

We define the time Th as the time when the kinetic energy per particle becomes comparable to the BEG transition 
temperature (so that the BEG is destroyed). Accordingly, employing Eq.(pll), 



g(0) ^ 



(24) 



Taking the ratio, we find 



.-1 

irr 
.-1 
h 



mT.Ll 



(25) 



where we have introduced the thermal length — \/rrCTc at the transition temperature. 

Thus, the loss of the phase memory can occur faster than the destruction of the condensate, only if the correlation 
length Lp is longer than ic [Q- From general physical considerations, it is clear that this condition is not realized 
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in the actual traps, if the main source of the noise is the background gas. The correlation length Lp is simply the 
thermal length of the hot gas, which is obviously much shorter than the thermal length of the confined cloud at Tc- 
It is also worth noting that similar conclusion can be reached in the general case, when the random field produces 
exchange between the components. Rigorous proof of this is also based on the Furutzu-Novikov theorem. 

The losses from the trap do not change the situation. The rate of the decay of the phase memory is practically 
comparable to the rate of the escape from the trap. Physically, this is so because the losses do not introduce any long 
range correlations which may disrupt the global relative phase faster than the loss of particles. This conclusion is 
based on the analysis of the corresponding Wigner equations following from the OPDM, where the collision integral 
is derived from the term (|l^). Detailed account of this will be presented elsewhere. 

It should be noted that a special attention should be paid to the situation of the quasi-equilibrium, when, e.g., 
the deposition to the trap and the evaporative cooling occur simultaneously, so that the trap remains at constant 
temperature and retains a fixed on average number of atoms. In this case, both rates of the deposition and of the 
losses can be large. Accordingly, the decoherence will be fast despite that no significant changes in the temperature 
and in the population are observed. The analysis of this situation will be presented elsewhere. Thus, if all the 
individual rates of disrupting the equilibrium are kept low on the experimental time scale, the correlations between the 
BEC components will persist on this scale. 

Below we will analyze the situation when the white noise is long ranged. As it will be seen, the above conclusion 
does not hold in this case any more. 



B. The phase decoherence induced by noisy external long ranged potential 

The noisy potential f could be induced by some external macroscopic fields. For example, some magnetic or 
electric fields can be imposed by extended macroscopic wires or by lasers. These fields fluctuate due to thermal 
and quantum effects. The magnetic trapping potential itself can fluctuate due to the fluctuations of, e.g., electric 
currents in the magnetic coils. The effect of the random term on the BEC dynamics has recently been analyzed in, 
however, different context in Ref. In such situations, the noisy field can be taken as ^(x, i) — U'{x.)ri{t), where 
f7'(x) stands for some macroscopic potential which breaks the symmetry, and ri{t) denotes the white noise factor 
obeying the relation 

{r,it)r,{t'))=AS{t-t'), (26) 

where the spectral weight A depends on the nature of the noise, and will be evaluated later in the case of the 
thermal equilibrium fluctuations. The Heisenberg equations for the two-component BEC become (we again ignore 
the self-interaction) 

ih^i.2 = [ ^ + C/±C/'77]*i,2, (27) 

where the sign is chosen for the first component, and "-" - for the second; U stands for the "symmetric" part 
of the trapping potential. 

Employing Eqs.(|2^), (p7|), and applying the Furutzu-Novikov theorem (see Appendix A ), we find the equations for 
the OPDM 

inpu{^,x',t) - _ [/(x) + u{k')-Zj[U'{^) - ;7'(x')]'}p.,(x,x',i); (28) 

*npi2(x, x', t) = _ [/(x) + [/(x') - *|^[C/'(x) + t/'(x')]'}pi2(x, x', t). (29) 



We will employ the above equations for estimating the magnitude of the ratio of interest Ti^jTh akin to Eqs.(p3|)-(p5[). 
The loss of the correlations between the components is determined by the term ~ A in Eq.(^), so that the rate 
Tj^r can be estimated as 

r^rl « f (30) 

where the bar denotes the effective averaging over the volume of the BEC. 

The heating rate can be found from Eqs.(p8|). We calculate the rate of the increase of the total energy in the trap 

as 
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^ = E / rfx|,,=,[^^ + [/(x')]p,.(x,x',t) = 

i 

E./ dx|,,=,Av2,(t/'(x)-[/'(x'))V..(x,x',i)«7V^[W^. (31) 



Accordingly, wc find from Eq.(]3l|) t,^ ^ « ^^[VJJ']^. Then, the ratio r-^l/Tf^ ^ becomes 



irr 
_-l 



(32) 



If t he pote ntial is long ranged as, e.g., U' ~ with some <5 > 0, and the cloud has a typical size Tc, then, 
C/'2/[VL/']2 w 7-^. Hence, we arrive at the estimate Eq.(|2^), where the role of the noise correlation length Lp is 
played by the condensate typical size Vc- Now it can be seen that, under the given circumstances, T^l/Tf^^ ^ 1, 
because Lc re- In other words, fluctuations of the trapping potential (or any other extended potential) destroy 
the inter-component coherence much faster than producing any substantial heating effect. This, however, does not 
mean that the rate of the decoherence is a significant issue in the actual traps, if the noise of the trapping potential 
is kept on the level of the thermal noise even at room temperature. Indeed, let us estimate the absolute value of the 
decoherence rate (|30|). For this purpose, we assume that the fluctuations ijU' are of the same order as those of the 
trapping potential U caused by the thermal noise 51^'^^ of the electric currents I^'^-' which create the magnetic 
trap. Then, we set 77 « W^^^ //(^), and impose the condition of the white noise: (W^'^) (t)(5/(^) (t')) = 2T/R5{t - t'), 
following from the classical limit of the fluctuation-dissipation theorem. He re R stands for the coil resistance. Then 
we find A — 2T/RI^'^^^ in Eq.(p^). For the mean in Eq.(pO|) we employ U'^ ^ ~ fi^, where ^iq stands for the 
mean chemical potential in the trap. Thus, Eq.(|30|) yields 

^ ^^7(^' ^^^^ 

where we have reinstated the Boltzmann constant ks- Eq.(p3|) gives — 10^^^ — 10^ ^'^s^^, for typical values 

T = 300K, fi/U — 10"^ — lO'^s"^, and for I^R — IW. Thus, the instrumental thermal noise does not produce any 
significant decoherence. 



C. Intrinsic decoherence of the global relative phase of the two-component BEC 

In this section we will discuss mechanisms of the decoherence induced by self-interaction in the trap. In principle, 
the two-component condensate can be viewed as a system in which certain information is imprinted into the relative 
phase between the components. Any quantum phase memory device, which could be employed in quantum information 
processing, should be able to retain the phase information for long enough time. Thus, it is important to understand 
the fundamental mechanisms of the relative phase memory loss. 

A well studied mechanism of the global phase decay at T — Q - the phase diffusion |^ can be understood in 
terms of fluctuations of the chemical potential induced by binary interactions and by a finite variance of the number 
of particles in the quantum coherent state. These fluctuations produce collapse of the global phase on some time Tc- 
It has also been found that the PD exhibits spontaneous revivals Q on some revival time tr N ^ Tc- The time 
ri{, which is too long on any experimental time scale, can be considered as a duration of the Poincare cycle. The 
spontaneous revivals indicate that the dynamics is reversible in a sense that it repeats itself. This, however, does 
not necessarily mean that the dynamics of 3> 1 particles can be physically reversed in time by imposing some 
physical external action. Thus, it is important to understand the physical reversibility of the quantum dynamics of 
the condensate on much shorter time scale Tc <t ^ tb,- This question will be addressed later. We will show that the 
time Tc is not the time after which the phase memory is lost irreversibly. The phase information can be recovered 
(partially) by imposing the time reversal pulse. 

Questions we are addressing here are: Can the evolution of the global relative be viewed as the PD process at 
T 7^ ? If it is the PD, is this process physically reversible? 

Before we start, it is worth noting that two conceptually very different situations should be clearly distinguished - i) 
the dynamics of the relative phase in the presence of the permanent exchange of bosons between the components, and 
ii) the phase dynamics with no exchange. In the case i), the relative global phase is characterized by some particular 
equilibrium value. Deviations from it increase the energy, and thereby activate irreversible processes of relaxation 
to the equilibrium. On the contrary, in the case ii) , if once prepared due to a short exchange of bosons, the global 
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relative phase cannot be viewed as a normal mode, if no exchange of bosons exists during the following evolution ( 
akin to the situation realized by the JILA group Q).We consider the case ii) in this paper. 

Here we will extend our analysis Ref. of the problem of the global phase time-correlations to the two-component 
case. The central element of this analysis is the concept of the projected Hamiltonian and the projected many body 
eigenstates, which will be often called just "projected states", of this Hamiltonian. The total numbers of bosons in 
the each component play a role of parameters (not the operators) in such a Hamiltonian, if no losses or exchange 
between the components exist^. 

We represent an exact many-body eigenstate |m, N, M), which is characterized by the following quantum numbers 
- the total number of bosons N — N2 + Ni and the half of the population difference M — (N2 — Ni)/2 as well as 
by a set of the quantum numbers m referring to the normal excitations, as an expansion 

\m,N,M)= J2 Cn^i,n^2....i^,N,M)\nio)\n2o)\nii)\n2i)... , (34) 

"-qO ^ Na~ {Uai + na2 + •■•): 

in the Fock space |?^lo)|^^2o)|f^ll)|n2l)... of the population numbers Uaj of some set of single particle states. Here the 
Greek index refers to the component a — 1,2, and the Latin index labels the single particle states in the corresponding 
component, so that Uao stands for the population of the ath condensate component. 

The expansion coefficients C„^j,„^2_...(r7i, A^, Af) form the Fock representation of the eigenstates. We call these 

coefficients the projected states, and introduce a short notation \m,N,M) for them pO| . 

It is important to note that, while being orthogonal with respect to the set m for given N, M, the projected states 
are noi orthogonal for different N, M. In other words, {M, N,m'\m, N, M) = 6m'm and {M' , N' ,m\m, N, M) 
Sm'M^N'n- Furthermore, in the case of the exact su{2) symmetry, the exact relation {M' , N, m\m' , N, M) = 
Sm'm holds regardless of the values M, M' . 

In 3D systems, close to equilibrium, the population factor of the condensate state in the each component is macro- 
scopic, and is characterized by relatively small fluctuations. If An^o stands for the fluctuation of the condensate 
population, and n^o is the mean population, the relative value Att-qo/'^qO ^ 1 p6[ |. This circumstance can be 
employed to express matrix elements of any condensate operator in terms of the overlap of the projected states. 
Speciflcally, let us say one wishes to find the matrix element (a2oaio)m,m' — {M,N,m\a\QaiQ\m' ,N,M — 1) of the 
condensate operators. Employing the representation ( ^^ and the smallness of the fluctuations of the condensate state 
populations, one finds 

(a2oaio)m,m' = Vnio(?i20 + l)C*„i,„„2,. ..(?", N, M)C„„i,„„2,...(™', N, M - 1) 

= y/^w^o[{M7N,m\m',N^I-l)+o(^^]]. (35) 

V "aO / 

Thus, the matrix elements of the condensate operators are given by the overlap of the corresponding projected eigen- 
states, with different values of the population difference. This simplifies considerably calculations of the correlators 
of the condensate operators. 

Indeed, the global relative phase time-correlation properties are completely described by the correlator Pi2{t) — 
(a|o(0'^2o(0'^2o('^)^io(0))- Employing Eq.(^, one obtains 

Pi2(t) = V"io"'2o(e ^ ' ^e ^ ' ') (36) 

within the chosen accuracy (here and below we employ units in which h ~ 1). H{N,M) stands for the pro- 
jected Hamiltonian defined as H{N,M) = J2m\'^^'^^'^)^rn{N, M){MJ^,m\ , with E^{N,M) being the ex- 
act many body eigenenergies ]lO| . The averaging is performed over the thermal ensemble. Explicitly, (...) = 
N M Prn{N, M){M, N , m\.. .\m, N, M), where Pm{N, M) stands for the canonical Boltzman factor. 



^The meaning of the projected Hamiltonian can be understood from the analogy with motion in a central potential, when 
the total angular momentum (and its z-component) is conserved. Due to the separation of the angular variables, the equation 
for the radial part can be formulated in terms of the Hamiltonian with some effective potential, which depends on the angular 
momentum as a parameter. This Hamiltonian for the radial part is the projected Hamiltonian. 
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The evolution operator g-^^f^-^^-i)* can be expressed as e-^-f^^^^^^^^^* = e-*^(^'^^)*Texp(-i di'iJ'(i')) , 
where the zeroth-order part is H{N,M) and the perturbation H' — H{N,M — 1) — H{N,AI) is taken in the 
interaction representation with respect to the zeroth-order part H'(t) — q>-H(n M)t j^i ^-iH(n M)t ^ Finally, Eq.(|3^) 
takes a form 



Pi2(t) = \^I^^(TexpH / dt'H'it'))), (37) 

Jo 

H{t) = e >{ __ + -^_^ + ...}e (38) 

within the chosen accuracy. It is clear that any correlator of the condensate operators can be represented in the form 
similar to Eq.(|^, (|3|). 

In fact, such an expansion is valid as long as a macroscopically populated BEC is present. Otherwise, the projected 
Hamiltonian is a very sharp function of M, because a transfer of even a single atom between the thermal components 
inevitably corresponds to creating excitations in the case of the broken su{2) symmetry. In the presence of the BEC, 
the each derivative term in the expansion of H'{t) introduces an additional factor 1/A^. Accordingly, it is enough to 
keep the first term only, so that H' — ^^^g^-^- H' can be viewed as an operator of the relative chemical potential 
= E,n{N,M) — E,n{N,M — 1) w dEm{N, M)/dM in the given eigenstates. The diagonal matrix elements of 
H' give the exact values of the The ensemble mean of it, then, becomes JI — '^^^ jy P{m, N, M)fijn = 

{dH{N, M)/dM) . It is important that H' = 0, if the intrinsic symmetry holds. Thus, another interpretation, 
H' is proportional to the change of the interaction energy per one transferred particle. As discussed in Ref. the 
proportionality factor is of the order of unity in the one-component BEC. In the two-component situation, this factor 
should contain the combination of the scattering lengths which vanishes at the point of the intrinsic su{2) symmetry. 

We employ the cumulant expansion of (M, N, m|Texp(— z J* dt' H' {t'))\m, N, M). Then, within the accuracy 1/N , 
we reproduce Beliaev's result Pi2{t) ^ exp{ijlt). We are, however, interested in the next non-vanishing effect 
of large but finite N. Therefore, the next cumulant should be retained as well, and the higher ones can safely be 
neglected within the 1 /iV^ accuracy. Accordingly, we write 

(M^,m|Texp(-i / dt'H'{t'))\m,'N^M) = 



Jo 

exp{i/i„t- /" dt' [ dT{M'^,m\ S H'{T)SH'{0)\m,'N^M) +o{l/N^)}, (39) 
Jo Jo 

where we have introduced the notation 5H'{t) for the off-diagonal part of the total matrix H'{t). The operator 
6H'{t) can be viewed as fluctuation of H' . 

As it has been discussed in Ref. ||l^], the double integral in the exponent determines the processes when adding or 
removing one boson to or from the BEC disturbs the normal component. In principle, a situation can be contemplated 
when such a disturbance is so strong that this term dominates, and rapidly suppresses the overlap of the projected 
states (differing in the total number of bosons by 1). This situation is similar to the orthogonality catastrophe (OC) 
occurring in Fermi liquids p8[ |. In Beliaev's work p7| , however, it has been proven that such processes are insignificant 
at T — 0. Thus, no the OC is to be anticipated in the BEC at T = 0. This implies that the double integral in 
Eq.(p9[) is finite and small in the limit t co. 

At finite T, without a rigorous proof, the Beliaev's result is widely applied to higher eigenstates jl^. Below, we 
will employ Eq.(|39|) and give semi-quantative general arguments which support the Beliaev's result at finite T. First, 
we note that the correlator in Eq.(^9|) can be calculated within the perturbation theory with respect to interaction 
between the quasiparticles. This has been done in the case of the single-component BEC in Ref. [0. The result of 
such calculations can be represented as 

^ F(m,7V,M)(M,iV,m|Texp{-i / dt'iJ'(t')}|m, TV, M) = e''^*"(*/^'')'-*/*oc , (40) 

m,N,M 







where 



r,-2 = l((M-7l)^>, toc«/ dt{5H'[t)5H'm. (41) 



The gaussian factor in Eq.(^) is a result of applying the central limit theorem for the averaging of e*'^"* over the 
ensemble. The quantity toe has been called the OC time [|l^. It determines how rapidly the overlap decays. Its 
decay is controlled by excitations. The rate f^J^ can be estimated as 
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t^l, « (5ff'2)r„ (42) 

where r„ stands for a typical relaxation time of the correlator {SH'{t)SH'{0)) which is taken in the exponential form 
^ exp(— t/r„). In 3D system the time r„ is finite because of very small statistical weight of the low energy modes. 
This time is given by a typical relaxation time of the normal modes characterized by the energies around the chemical 
potential. Obviously, r„ <^ Td (|4l]). Thus, taking into account the nature of the operator H' discussed above (see 
below Eq.(^)), we conclude that ~ 1/-^, because {SH' {0)SH' (0)) ~ 1/N as fluctuation of any extensive quantity 
taken per one particle. It should, however, be noted that the form ( ^0| ) is not valid in the limit t ^ oo. In Rcf. [ p^ , 
we have discussed that the double integral in Eq.(|39|) should remain constant and small in the limit t = oo {t> toe) 
even at T 7^ 0. The reason for this is the following: the spectral weight of the correlator {SH' {t)6H' {0}) in the limit 
of zero frequencies is effectively collected from the lowest order processes of scattering of the quasiparticles with almost 
zero energy transfer. The term H' , which is actually proportional to the interaction energy, removes the degeneracy 
between the corresponding quasiparticle states. This suppresses the spectral weight in a very similar manner to the 
effect of the level repulsion studied in the Random Matrix Theory This is essentially a non-perturbative effect, 
and it becomes important at frequencies smaller than a typical value of the interaction matrix element, which is scaled 
as ~ 1/VN Thus, the exponential factor q-^I^oc \^ Eq.(^0|) is correct only for times shorter than some time 

t' ^ \/iV, and it levels off to a constant ^ exp(— i'/ioc) ~ l + o(l/\/]V) (note that toe ^ N ) as t becomes longer 
than t' . 

The quantity Td is determined by the ensemble fiuctuations of the chemical potential, and it describes the rate of 
the phase diffusion at finite temperature. We will call it the thermal PD in order to distinguish it from the quantum 
PD 0. Thus, tT^ - l/^/]V > ~ and, therefore, the OC is irrelevant because the reversible dephasing 

dominates in Eq.([lO|) on any practical time scale, provided no extrinsic factors play any role. 

The evolution of the condensate operators can be viewed as a result of the ensemble averaging of the non- decaying 
exponents, each representing a particular realization of the chemical potential. This is so, if the relation 

{A'Cn, m\m', N^I ± 1) = (5„,,,„, (43) 

holds. We note that Eq.(|4^) is, in fact, an extension of the Beliaev's result for finite T formulated in terms 
of the overlap of the projected states. Accordingly, the correlator (^6|) simply becomes 

Pi2{t) = ^/^^Q P{m,N,M)e'^'-' ^ V^wn^e'^'^^'/^^^" (44) 

m,N,M 

within the central limit approximation. 

The dephasing time Td can be calculated in any order with respect to the interparticle interaction. For this 
purpose, in the weakly interacting system, it is enough to consider the Bogolubov gas of the non-interacting quasi- 
particles as a model of the excited states. Specifically, the eigenenergies are represented as Em{N, M) = Eo{N, M) + 
J2a^a{N, M)na, where Eo{N,M) is the ground state energy and ea{N,M) stands for the spectrum of the quasi- 
particles characterized by the population factors na- Then, Eq. (p|) yields 

rd' 

where Tc stands for the quantum phase diffusion (at T — 0) collapse time. The quantum PD has been analyzed 
by many authors in Refs. [|6| in the case when the initial state is the coherent state with well defined global phase. 
We note that, in the two-component BEG, the shot noise in the initial deposition of N should also contribute to 
the dephasing of the condensate correlators at T = even in the situation when no initial global relative phase 
existed. Furthermore, in actuality, the shot noise is most likely to wash out completely the effect of the quantum PD 
in Eq. (^6|) . In our work ||], we discussed this in detail, and have found that in *^Rb the collapse time Tc can be 
as short as 30ms, if the shot noise reaches values about 10%. This should be compared with typical estimates of the 
quantum PD collapse times ~ 1 — 10s 



* that the action of the condensate operators on the ground state does not produce excitations within 1/7V accuracy 
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It should also be noted that = in the case of the exact su{2) symmetry. Here we will not focus on calculating 
the dephasing rates for the specific geometries and will use it as a parameter. 

Summarizing, the evolution of the condensate operators does not disturb the normal component, and, as a conse- 
quence, the correlators of the condensate operators exhibit reversible dephasing determined by the ensemble fluctua- 
tions of the relative chemical potential. In the next section we will discuss how such a reversibility can be tested in 
the atomic echo experiment. 



IV. ATOMIC ECHO 



As discussed above, the evolution of the condensate operators in 3D system including finite number of bosons 
does not disturb significantly the normal component. The decoherence of the condensate operators is a result of the 
ensemble averaging, with no decoherence occurring in any eigenstate. In this section, we will address the question of 
physical reversibility, so that an initial phase information can be recovered on times longer than the dephasing time 

Td- 

In general, the echo effect (spin echo, photon echo ||2^], plasma echo 22 1) - stimulated revival of some quantity ~ is 
widely employed as a test for reversibility of dynamics. Once initially prepared, the quantity may exhibit an apparent 
decay, which is a sort of inhomogeneous broadening without actual loss of the memory of the initial state. Then, if 
some time-reversal pulse is imposed at later time T , which can be much longer than the dephasing time , the 
quantity may revive completely or partially at time t ~ 2T, and, if the time reversal pulse is strong, at integers of 
T (multiple echoes). 

Generally speaking, any Hamiltonian dynamics is reversible in time. This means that the formal time-reversal of 
the evolution will result in the restoration of the initial state. This, however, cannot always be done physically, that 
is, by applying some physical disturbance to the system. Thus, a proper criterion for the physical reversibility of the 
dynamics should be introduced. We will employ the criterion which relies on the strength of the echo. Specifically: 

// an external pulse can induce the echo whose strength is significantly larger than fluctuations (quantum or statis- 
tical) of the measured quantity, the evolution of this quantity can be considered (partially) reversible. 

It is important to realize that the echo in the condensate can be observed even though the time delay T is much 
longer than typical relaxation times of the excitations in the system and/or in the bath. The time delay should be 
shorter than the irreversibility time Ti^r (^3|) , ( |3C| ) . As discussed in Sec. Ill, this time is determined by the strength 
of the extrinsic fluctuations, which can even be treated within the (5-correlated approximation, that is, within the 
assumption of infinitely short relaxation times. 

In our work we have presented the echo solution in the two-component BEG in the case of the shot noise 
dominated dephasing. In this situation, the echo strength reaches 100%, and is achieved by imposing the 7r-type time 
reversal pulse. Here we will analyze the echo in the case of the intrinsic dephasing, that is, when either the quantum 
or thermal PD takes place. 



A. Atomic echo at T = 

Let us consider the echo experiment in the context of the three pulses scheme We take Ui as the evolution 
operator which corresponds to the tt/2 pulse which converts the initial state with N2 ^ N, Ni — into a coherent 
state of the two condensates Then, some time-reversal pulse U' is imposed at time T > later, and, finally, 
the read out 7r/2 pulse ||] is imposed at time t > T . 

We choose, as a measured quantity, the population difference represented by the operator Iz in eq.(|^) and 
measured by the JILA group [^. A state vector is now labeled as \j, M) by the conserving total angular momentum 
j = iV/2 and by its z-projection M . At T = 0, the excitation label "m" can be suppressed. Then, the mean after 
the read out pulse becomes 

(4) = (j,j|C/ie'™([/')te'(*-^)^C/tj,[/ie-(*-^)^C/'e-™[/|b-,j) . (47) 

The pulse operators Ui, U' can be found in the sudden approximation, if the duration of the pulses is short. Then, 
a general pulse operator is ?7(|f |, 0) = cxp(— ?(«*/+ + iJ.c.)), where v+ — J dtfl(t) — \v\ exp(i(/)) is the time integral 
(over the pulse duration) of the Rabi frequency Q{t) in Eq.(|l|), and the operators I± are represented in Eq.(||). The 
value \v\ is the magnitude of the pulse, and (f> stands for its phase. Specifically, for the 7r/2-pulse, one can choose 
Ui — U{n/A, (j) = — 7r/2) = exp(i^/y). For the time reversal pulse, we will take U' = exp{ipiy), where (3 stands for 
the pulse strength. The dependence of the echo strength on (3 will be investigated below. 
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It is worth noting that the operator e*''^" is the operator of finite rotations, and /3 stands for the corresponding 
Euiler angle. An exphcit form for the matrix elements d!]^), ^{(3) = {M' , j\e^^^y |j, M) is well known for any arbitrary 
dimension of the su{2) representation (see, e.g., in Ref. p3[). For example. 



4^.W = a.|e-1.,M) - (eos0"" (sm0^"'' (48) 

where M = —j, — 7 + 1, — The expression for arbitrary M,M' involves Jacobi polynomials, and can be 
found in Ref. |2^ . 

Thus, the mean ( ^ ) can be calculated as a finite product of the matrices + 1 by N + I. In our work |^ we 
have analyzed the exact echo solution, when the main source of the dephasing is the shot noise. In what follows, we 
study the echo effect in the situation of the quantum PD As it can be seen later, the echo solution can be found, 
practically, exactly as well. 

First, we remind that, under the given initial condition, the evolution of the population difference represented by 
Eg. (^7]) exhibits collapses and spontaneous revivals [||. Indeed, setting the time reversal pulse to zero in Eq.(^7|), one 
finds 

{Iz) = {j,M^'^e^'"j.e-''"e-^^'^\j,j), (49) 

where we have employed the standard su{2) relation e^^^-^H/^e'?'^'' ~ I^- Taking into account Eqs.(^), and the 
dependence of the chemical potential on the population difference, one immediately obtains the collapses revivals 
solution where the collapse rate is ( ^6|) is a factor of \/n larger than the revival rate . In what 

follows, we will assume that A^ ^ 1, so that the revival time is too long on any practical time scale. Furthermore, 
as we mentioned above, even weak shot noise should wash out completely any spontaneous revivals. Thus, we will be 
always assuming that Tc < T <^ t-ji, so that the echo - stimulated revivals would not be confused with the spontaneous 
revivals. 

Calculations of the mean ( ^7|) exactly for large j and for (3^0 is a very combersome problem. We simplify it by 
noting that, as indicated by Eq.(^), the initial 9 = 7r/2-pulse creates states with |M| « ^/J <^ j. Accordingly, the 
time reversal pulse U' acts effectively between states with such M . This, implies that the formalism of the quantum 
relative phase can be applied for calculating the matrix elements of e^^^^ . We introduce the phase 1^9 as a conjugate 
variable to Iz- Thus, the representation of the operators (^ can be chosen as 

g 1 d d 

Iz = -i-g^ , Ix ^ jcosip + ——cos(p— + o{l/f), 

(50) 

J__5_ . d_ 

2j dip ^ dtp 

which is a formal expansion in 1/j = 2/N. The terms shown explicitly satisfy the standard su{2) commutation 
relations within the accuracy 1/N'^. We emphasize that the quantum phase representation is valid only on the space 
of the states characterized by |Af| < -^/J. As it will be seen later, the maximum echo can be achieved for the strengths 
of the pulses which satisfy the relation (3 ~ l/v^- Accordingly, within the formal series (|50|), (3Iy « pjsmip and 
^ipiy Qi(3jsmip^ gQ ^-f^^ matrix elements become 

d^kM « Jm-m' m, \M\, \M'\ « J » 1, (51) 

where J„(a;) is the Bessel function of the first kind of the nth order. This relation can also be obtained from the 
hypergeometric series, which yields both the limiting representation of the Bessel functions, and the Jacobi polynomials 
p^ . This approximation simplifies calculations of the mean ( ^7| ) considerably. On Fig.l, we have compared the results 
of the exact numerical calculations of (|4^) with the calculations within the quantum phase method. As can be seen, 
agreement is very good (for N w 100), and it is naturally expected to become better as A^ increases. 

Accordingly, in what follows we will be discussing the quantum phase method only. 

Thus, the complete representation of the mean ( |47| ) becomes 

(/.)=2Re J2 d%{^)d^ll^{-^)iI,)KK-lJK-MmJK-M'-im^ 
M,M',K 

g il{Eo(M)-Eo{M'))T+(E„{K)-Eo(K-l)){t-T)]^ (^^2) 

Eo{M) = (eo + hN)M + ^M^ + ^N^, (53) 



ly = j sin p + — — sin p— h o(l/j" 
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where {Ix)kk-i = {Ix)k-ik = | + K){j — + 1) stands for the matrix elements of the operator I^; we have 
introduced the ground state eigenvalue Eq (M) of the Hamiltonian (|^) as a function of the half of the population 
difference M and of the total number of bosons TV, with 61.2,3 being the coefficients which can be obtained 
from Eq.(0) within, e.g., the two- mode approximation [ p5| . The results of the numerical calculations of Eq.(|5^) are 
presented on Fig. 2 (for bi = 0). As can be seen, the multiple echoes occur at times t « kT, k — 2,3,4, .... 

The multiple echo solution can be found, practically, exactly in the limit T ^ Tc ~ 2(^/762)^^. As it will be 
seen later, the maximum echo is achieved when /3 « rc/(T-\/j), so that the argument of the Bessel functions is 
~ \/jTc/T ^ y/j. The Bessel function Jn~m{z) becomes essentially zero when its order |n — m| > \z\ |2^. Thus, 
in Eg. (|5^) , while the index M runs over the range w [— a/J; the other indices stay relatively close to M. 

Accordingly, in the sum (^2|), we approximate (j + K) {j — K + 1) fa j. Then, the sum over K can formally 
be extended from —00 to +00. This, allows to calculate this sum in a closed form by employing Eq.(|5^) and by 
implementing the identity (9.1.79) of Ref. |Q. Accordingly, Eq.(^2|) acquires the form 

(/.) - J E JL+im,it - r))e-(*-^(^+^»Vr.= eos[eo(t - r(L + 1)) + ^^^^^], (54) 

where we have introduced the notation L = M' — M , and have eliminated all the terms with L < because their 
contributions become significant only at negative times. We also employed the central limit theorem and made the 
replacement Y.M d%ij)d^M+Lji-l) exp[i62(M + L/2){t-T{L + 1))] w exp[-(i - T{L + l)f/T^], which is valid for 

Eq.(Q) describes the stimulated revivals - the echoes. They occur at the times t « T{L + 1), L> 1. The analytical 
solution (|4|) is shown on Fig. 2 together with the numerical one, both calculated within the quantum phase method. 
The agreement is quite good. 

It is interesting to note that the magnitude of the echoes exhibits non-monotonous dependence on the strength 
f3 of the time-reversal pulse. Another feature is that the maximum echo can be created by the weak time-reversal 
pulse. Let us consider this in detail for the first echo (L = 1) at the maximum of the gaussian, which occurs at 
t = 2T. Then, the relative echo magnitude becomes 

— = J2{Pjb2T) = J2{2f3y^T/T,), (55) 
J 

where we have expressed 62 in terms of the collapse time, which follows from Eqs.(^), ( p3| ) as = \/j|52|/2 (we 
ignore the shot noise). Taking into account that the first maximum J2{z) ~ 0.5 occurs at z « 3 |2j], we find the 
magnitude of the strength of the time-reversal pulse 

required to achieve the maximum echo. Conversely, J2{z) — for z « 5 for the first time, and the echo vanishes for 
(3 « ^.iPmax- Taking the higher value of (3 results in the oscillatory dependence of the echo strength on (3, with the 
total amplitude slowly diminishing as ~ \l\f^. Specifically, for /3 ~ 1 , Eq.(^5|) yields \lz\l 3 ^ i^^^* <C 1. As it 
will be seen later, at T ^ this dependence can cross over to \lz\l 3 ^ J^^^^ ^ li which makes the echo practically 
zero because it is on the level of the statistical noise. 

The above result indicates that, in order to revive the quantum phase, which exhibits the quantum phase diffusion 
, it is enough to exchange coherently the number of bosons « /?j ^ \/N between the components. 
In the next section, we will consider the effect of finite temperature on the echo. 

B. Echo at T / 

As discussed in Sec.IIIC, the condensate operators acting on the eigenstates do not mix them with other eigenstates. 
This is stated by Eq.(B), which is an extension of the Beliaev's result jl^ obtained for T = 0, and which is widely 
employed at T 7^ |19[. Nevertheless, a special attention should be focused on calculating the matrix elements of 
the time reversal pulse e*''^" at finite T, when the normal component is present. Indeed, a formal separation of the 
operator ly ^ N into the condensate and the normal parts could be expected to produce a sort of the Debye- Waller 
thermal factor exp(— /3^A^'), where N' stands for the number of particles in the thermal cloud. This logic, however, 
does not take into account the symmetry considerations. 
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We start our analysis of this problem by noting that, in the case of the exact su{2) symmetry, the eigenstates 
|m, -/V, M)(o) of the symmetric Hamiltonian i/(o) form the su{2) representations for every value from the set of the 
excitation quantum numbers to. This set should contain also the value j characterizing the dimension of the corre- 
sponding representation. To emphasize this, we will use the notation \m, j, N, M) for the state |to, A^, M)(o) char- 
acterized by some value j. If in the initial state, all the bosons belong to one component (as in the JILA experiment 
1^), this value is j — N/2, and it remains unchanged during following evolution. This is so because the Hamiltonian 
forms a closed algebra with the generators (Q) . Accordingly, 

(M, A, J, m\e'0'y |m', j, A^, M') - d^^,, (/?)<5„,„, (57) 

for the initial set of states characterized by some j. 

The situation is different if the intrinsic symmetry is broken, so that the Hamiltonian does not form a closed algebra 
any more. Then, it is natural to expect that the value j is not a good quantum number under this circumstance. 
Furthermore, it is clear that, at T > Tc, this value should relax to j ^ VN, even though initially it was j — N/2. 
Indeed, the square of the total angular momentum operator (||) is = /+/- + — Iz- At the point of equal 
populations, the means (Iz) = 0, (J^) ~ A^ and (/+/_) A^, where we have employed Eq.(|^) and have taken into 
account that no long range order exists above T^. This should be contrasted with the situation at T = 0, when 
practically all bosons occupy the condensate states, and (J_|_/_) ^ N^. 

What happens at T ^ and T < Tc7 The dimension of the dominant representation is given by the 
numbers of the condensate atoms. Indeed, keeping in mind the explicit form (|^), one can write (/_|_/_) — 
(/dx^'^(x)5'i(x) /dx'*|(x')^'2(x')) = A^c2A^ci(l + o(l/A^)) « {N/2)'^ for equal populations in the uniform case. 
Furthermore, the correlator (/+(t)/_(i)/+(0)/_(0)) remains time independent as long as the BEC is present. Indeed, 
employing Eq.(|43[), it is easy to show that this correlator does not depend on time in the main ^ N"^ limit. Thus, 
the dimension of the representation is practically selected by the initial condition. 

Now let us consider how Eq.(^) changes when the su{2) symmetry is broken intrinsically, provided the population 
difference M is still conserved (due to the absence of the exchange between the components). We employ Eq.(|5~' 
in calculating the matrix element and find (M', A, TO'|e*''^« |to, A^, M) = J2n,j ^nWi^'^' ^ ^^^'l"^' 

{M,N,j,n\m,N,M). This is equal to J2nj '^M\MW{MO^,rn'\n, jJV, M') x {MJ^,j,n\m,'N^M) due to the 
property of the projected states that their overlaps are equal to the overlaps of the corresponding eigenstates for the 
same values of M, A^ |^. Finally, due to the orthonormality of the projected states in the subspace with given j, we 
find 

N/2 

{M',N,m'\e''^^y\m,N,M) = ^d^l^^j^j,{l3){M' ,N,m'\Vj\m,N,M), (58) 

3 

where Vj projects the eigenstates to the space of functions with angular momentum j, and the lower limit is either 



X 



or 1/2, depending on whether A^ is even or odd, respectively. Eq. (p8|) is a generalization of Eq.(57). This equation 
gives the representation of the matrix elements of the pulse operator in terms of the universal values d^l) (/3) |2^ ] 
and the overlaps of the exact projected eigenstates. Employing Eq. (|58|) as well as the orthonormality of the projected 
eigenstates for the same M, we find the following representation for Eq.(^) 

31,32,33 Ml, 2, 3, 4 

G{Mi,M2,M3,M4,t) (59) 

where j = N/2; {Ix)^m^ Ms stands for the matrix elements of the operator calculated in the space of the angular 
momentum j2] and 

G (Mi,M2,M3,M4,i) = 

Pn^{N){jXm\e'^"^"^^V,,e'^'-^^"^'''^^V,^^^^^ , (60) 



^ It is important to note that the projected eigenstates \m,j, N, M) do not depend on M. To continue the parallel between 
the projected states and the radial part of the wave function in a central potential (see Sec. Ill C), we note that this radial part 
does not depend on the ^-projection M of the angular momentum. 
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where we have chosen the initial state characterized by j = N/2 (initial populations N2 = N, A^i = |^) and have 
suppressed the parameter N in the notation H{N,M) of the projected Hamiltonian; Pm{N) denotes the canonical 
normalized thermal distribution of the initial states (characterized by M — j, j — N/2). 

We note that Eqs.(|59|), ( |60|) are exact. They, however, can conveniently be employed only if the BEC is present. As 
discussed above, in the case of the dominant population of the BEC state, the dimension 2j + 1 of the su{2) repre- 
sentation is dictated by the initial condition N2 = N, Ni — 0. Accordingly, one can choose ii = J2 = J3 = j = N/2 
(that is, omit the external summation in Eq.(p9|)), and set Vj — 1 in Eq.(|60|). At T = 0, Eq.(]59|) naturally transforms 
into Eq. (^) , if one chooses the only term m — , that is the ground state, in the sum (|60|). 

In what follows, we will only estimate the factor ( |60| ) for the first echo within the cumulant expansion, which was 
employed in Sec.IIIC and is outlined in Appendix B. We consider temperatures low enough, so that the fluctuations 
of the condensate populations can be ignoredj^. 

The suppression of the echo at t = 2T is described by the factor (|60|), where we set t = 2T. Employing Eq.(p9|) 
as shown in Appendix B, we find 

( h{2T))^j 4fMW2)4Vi,,(-V2)JM'-M(/3j)JM'-M-2(/3j) X 

M,M' 

Y^^^p^^(^j^,jyT{E„,(M~l)-E„,(M)+E^(M')^E^(M' -l)-[2(M-M' ^^g-^^ 

where we have taken into account that the effective values of M, M' are much less than j = N/2] the definition 
of the effective distribution function Pm(M) created by the first pulse is given in Appendix B, Eq. (p3|) . We have 



approximated the matrices of finite rotations by the Bessel functions similarly to how it was done in Eq.(52). The 
matrix has also been replaced by j /2 in this space, and, accordingly, the terms M3 = M2 ± 1 have been 
selected. 

We note that Eq.(|5^) (with t = 2T ) is, in fact, the limiting case toe — > 00 of Eq.(|6ll). Thus, we have actually 
assumed that the limit toe ^ 00 corresponds to T — 0, and, conversely, toe becomes small enough to include its 
effect at high T. Here we will not calculate a specific expression for toe as a function of T. 

In order to obtain a closed form expression for Eq. (|6l|) , we employ the smallness \AI/N\ <C 1, \M' /N\ <C 1, and, 
accordingly, expand the eigencnergics as 

EUM) = E„,{0) + e„.Af + ^62™^' + 0(1/ N^), (62) 

_ dE„,{M) _(PEUM) 

oM \M=a oM'^ \M=o 

around the point M = of equal populations. Then, Eq. (|6l|) becomes 

(/,(2r))«jX^^(/3jVK-2(/3j)e'^^^^^"['(^-^)'+^l*«^>, (64) 

K 

where we have ignored the fluctuations of 62m, ~ 1/-/V and wrote P^iM) exp[iT{M' — M)b2m] = exp[iT(M' — 
M)&2], with 62 = J2m-Pm{0)b2m- It is worth noting, that, in the hmit toe — 00, Eq.(|6^) yields the first 
echo magnitude ( |55| ) (where the value 62 = fe2m=o is replaced by 62). We have also taken into account that 
J2m '^j"'!/ ('''/2)'^m+i j(^/2) = 1 with a very good accuracy in the limit j ^ 1. 

The exponential factor in Eq.(64), containing tQ^, produces some suppression of the echo strength. As can be 



verified numerically, the sum ( |64| ) at its maximum decreases from the value about 0.45, which corresponds to toe = 00, 
to the value about 0.06, corresponding to toe ~ 2T. A change is also exhibited by the value of /3j at the maximum 
of the echo. It shifts from the value given by Eq.(^) (3 « (in the hmit toe = 00) to /3 ~ 4.5/j for toe ~ 2T. 

This can be understood as follows, as toe ~^ 2T, only few Bessel functions depending on the product (3j contribute 
to the sum (|4|), and /3j « 4.5 corresponds to the maximum of this contribution. As we discussed in Sec.IIIC, 
toe ~ -/V as long as the population of the condensate is comparable to N. Thus, practically, toe ^ T Ki Td ^ ^/N. 
This implies that the thermal effects do not suppress the echo substantially, if compared with the zero T limit. It 
should, however, be noted that at temperatures high enough, the fluctuations of the condensate population [ p^ would 
lead to the necessity of including the summation over ji,2.3 in Eq. (|59|) . This problem will be analyzed elsewhere. 



Otherwise, the complete form ( [59| ) with the summation over ji,2,3 should be analyzed 
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We also not that external factors, which introduce irreversible decoherence, may suppress the echo significantly. If 
a typical time of the loss of the coherence due to these factors is Tirr, as given by, e.g., Eq.(^, the first echo intensity 
will acquire the exponential factor ~ g-zr/xi,.^^ ^^^^ accordingly, no echo will be practically observed in the situation 



V. SUMMARY 



We have shown that the intrinsic symmetry in the multi-component atomic mixtures has critical impact on the 
decoherence of the Ramsey fringes. This symmetry should be broken in one way or another in order to induce 
damping. The extrinsic factors, such as, e.g., the interaction with the background gas, can produce the decoherence 
of the inter-component correlator. This decoherence, however, is not faster than the rates of the heating and of the 
induced losses. 

We have shown that, in the two-component BEC, the intrinsic decoherence is a result of the ensemble fluctuations 
of the relative chemical potential. In other words, the evolution of the condensate operators occurs as though the 
normal component is not affected by it. Therefore, this evolution can be considered reversible. 

We have analyzed the atomic echo effect in the Ramsey spectroscopy of the two-component BEC. The strength of 
the echo is obtained in the case of the phase diffusion at zero and finite temperatures. We have found that the echo 
survives finite temperatures, as long as the population of the condensate state remains dominant. 

We also point out that the echo effect depends essentially on the nature of the many-body correlations, and as 
well as on the deviations from the intrinsic su(2) symmetry. Accordingly, its measurement would provide valuable 
information on the nature of the many body correlations in the trapped multi component BEC. Because of this, we 
suggest that the atomic echo effect should be studied experimentally. 
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APPENDIX A: EQUATION FOR OPDM IN THE PRESENCE OF WHITE NOISE 

Furutzu-Novikov llial theorem states that 



{ax,t)zix\t')) 



dy / dr(e(x,t)e(2/,r))( 



(Al) 



where the averaging is performed over the gaussian noise £^{x,t); Z{x,t) is an arbitrary smooth functional of ^. 
This relation can be verified by expanding Z(x,t) in the functional powers of ^ and by comparing the l.h.s. and 
the r.h.s. ter m b y term. 

We apply (Al) to the Heisenberg equations, following from Eq.(^ and written for the operators py(x, x',i) = 
M/t(x,t)M/,(x',t). 



inp^Jx,x',t) = [-iJox - 6(x,t) + Fox' +Cj(x',t)]pjj(x,x',i), ^1 = -^2 = 



(A2) 



where no two-body interaction is taken into account. We average Eq.(A2) over the initial condition and the noise, 
and obtain for the OPDM /9y (x, x',i) — (p^ (x, x', t)): 



ihpij{:x.,x.',t) = [-iJox + i?ox']/3y(x,x',i) -I- ([-^i(x,t) + ^j(x,t)]p,y (x,x',t)). 



(A3) 



Then, we represent the mean ([— ^i(x,i) 4- ^j(x, t)]pij(x, x',t)) in accordance with Eq.(Al), and write down the 
equation for the functional derivative as 



ih 



' = [-Ho. + gox'] : + ([-^^(x, t) + ej(x, t)]- ^ 



d (5pt.,-(x, x^t) 

+ [-<5(x - y) + ^(x' - y)](p., (x, x', t))S{t - r) 



(A4) 
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We take into account the causality, so that 5pij{x.,^x' ,t)/6^{y,T) = for t < t. Then, we find 



6pij{:x.,x',t) 



from Eg . (|A4| ) . We also take advan tag e of the white noise structure (|l|). Finally, substituting Eq.(|l|) and Eq.(^) 
into Eq.(|Al|), we obtain Eqs.(|6|), (|l|). 

APPENDIX B: THERMAL SUPPRESSION OF THE ECHO 

To simplify the notations, we will be suppressing N - the total number of bosons. We take ji,2,3 = j, and represent 
the factor (^0|) at i = 2T , and for M3 = M2 — 1 (the term with A/3 — M2 + 1 is just the complex conjugate) as 

G(Mi,M2, Af2 - l,M4,2r) - ^P„-(j;^'|m^4)(M4r?^l"Aj') x 

and select the diagonal terms (with respect to the excitation label) only m — n. In fact, the evolution of the off- 
diagonal terms is controlled by the excitations, and, therefore, they decay rapidly on the chosen time scale t^. Then, 



Eq.(Bl) acquires the form: 

G{Mi,M2, M2 - 1, M4, 2T) = 

m 

,m\m',j). (B3) 

m' 

The quantity PmiMi) can be considered as an effective distribution function created by the first pulse for the given 
population difference 2A/4. As it can be seen, it satisfies the normalization condition X]m^"i(-^^4) ^ ^■ 

As seen from Eqs.js^), (^), the first echo is dominated by the term Af4 — Mi + 1. Accordingly, in what follows 
we will consider the term G{M — 1, AT', M' — 1, Af, 2T) only. We apply the cumulant expansion in order to evaluate 
it. Thus, 

G{M - 1, A/', M' - 1, A/, 2T) - ^ P^(A/)(Af^|e-*^^(^^)e'^^(*^-i) x 

m 

e ^™''^e-'^"^^''-'^\n^I) = Gi,G2,. (B4) 

The first cumulant 

Gic = P,„(Af)e^^[^'"(*^-i)-^™(*^)+^™(*^')-^'"(*^'-i)l (541) (B5) 



corresponds to the assumption (AT, m|m', Af) = Sm,m' for any value of the difference A/ — A/' akin to Eq.(^3|). 
This, however, is not completely the case. In fact, for large values of AT — M', a sort of the Debye- Waller factor 
enters this product. The second cumulant G2C takes care of this effect. 
We employ the identity 



g-zff(M')T ^ e~^^(^'^)^Texphi(Af' - M) dt'H'{t% (B6) 



^Texp[-i(Af' - Af) [ dt'H'{t')] 
Jo 



H'{t) e'«(*^)*{^f^ + o(l/^)}e-*^(^'^'*, (B7) 

where we have followed the same reasoning which lead us to Eq.(^) - keeping the first term of the expansion of the 
projected Hamiltonian. 

We are interested in evaluating the effect of suppression of the diagonal matrix elements (with different M and 
the same m) on the time evolution. Then, we write (Af, m|e"*^^(^^)e''^^(*^^^)e*'^^(^^')e^*'^^(^^'"i) |m, Af) « 
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(M,TO|e'™(^^-i)|r7^7M)(Af^|e-'™(*^)|r?CM)(i\^ which after employ- 

ing Eqs.(|B(]|),(B7) for the each mean, and substituting back to Eq.(^J), yields the first (B5) and second cumulant 

G2c = exp{-[2(M-M'-i)2 + ^ ] x 



dt' / dtY.^P„,{M){M,N,m\SH'{t)SH\0)\m,N,M)}, 



(B8) 

where SH'{t) is defined below Eq. (|3^) . We employ the definition of toe in Eq. (^) , and rewrite Eq.(B4) as 

G (M - 1, M', M' - 1, A/, 2T) = 

J2 F,„(M)e*'^<^"(*^~i)~-^'"(^^)+-^"(^^')~-^"(^^'"i)'"[2(M-M'-i)2+|]t^J,}^ j-gg^l 

which is valid in the limit Tn <T <^ toe ■ This expression is to be employed in Eq. (B9h . 
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Fig. 1 
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FIG. 1. Comparison of the results of the exact calculations of Eq.(|47|) with that obtained within the quantum phase method. 
Solid line corresponds to the quantum phase method. Dotted line - the exact calculations. Parameters chosen: j = 45 (total 
number of particles AT = 2j = 90), /3 = 0.22, T = 10, 62 = 0.05, eo = 2. 
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FIG. 2. The echo effect in a confined two-component Bose-Einstein condensate at T — 0. Vertical axis is labeled by the 
relative magnitude of the population difference after the read-out pulse: Iz is the half of the population difference; j stands for 
the half of the total number of bosons. Solid line - the analytical solution (^); Dotted line - the numerical solution of Eq.(^2[). 
Parameters chosen: j = 2000 (total number of particles N = 2j = 4000), T = 2, 62 = 0.1, eo = 40; /? = 0.0005, f3j = 1 
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Fig. 2b 
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